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Abstract. The spectrum of a quantum system has in general bound, scattering and resonant parts. 
The Hilbert space includes only the bound and scattering spectra, and discards the resonances. One 
must therefore enlarge the Hilbert space to a rigged Hilbert space, within which the physical bound, 
scattering and resonance spectra are included on the same footing. In these lectures, I will explain 
how this is done. 
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INTRODUCTION: LECTURE 1 

t> 

In Quantum Mechanics, observable quantities are represented by linear operators. The 
eigenvalues of an operator represent the possible values of the measurement of the corre- 
sponding observable. These eigenvalues, which mathematically correspond to the spec- 
trum of the operator, can be discrete (as the energies of a particle in a box), continuous 
(as the energies of a free, unconstrained particle), resonant (as in a decay), or a combi- 
nation thereof. 

The Hilbert space includes only the bound and scattering spectra, because the Hilbert 
q-i space spectrum of an observable is real, thereby discarding the resonance spectrum as 

unphysical. However, radioactive nuclei and unstable elementary particles are physical 
objects that ought to have a place in the quantum mechanical formalism. This is why we 
need to extend the Hilbert space to a rigged Hilbert space, within which the resonance 
spectrum has a place. 

The purpose of this series of lectures is to explain how one should use the rigged 
Hilbert space in quantum mechanics and, in particular, how to incorporate the resonance 
spectrum into the quantum mechanical formalism by using the rigged Hilbert space. 

When the spectrum of an observable A is discrete and A is bounded, then A is defined 
on the whole of the Hilbert space Jrf? and the eigenvectors of A belong to J^. In this 
case, A can be essentially seen as a matrix. This means that, as far as discrete spectrum 
is concerned, there is no need to extend M '. However, quantum mechanical observables 
are in general unbounded and their spectrum has in general a continuous part. In order to 
deal with continuous spectrum, we use Dirac's bra-ket formalism. This formalism does 
not fit within the Hilbert space alone, but within the rigged Hilbert space. 

Loosely speaking, a rigged Hilbert space (also called a Gelfand triplet) is a triad of 
spaces 

$CJfC<I) x (1) 
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such that ffl is a Hilbert space, 4> is a dense subspace of J^, and 4> x is the space of 
antilinear functionals over 4>. Mathematically, 4> is the space of test functions, and 4> x 
is the space of distributions. The space 3> x is called the antidual space of 4>. Associated 
with the rigged Hilbert space (1), there is always another rigged Hilbert space, 

•DcifC*', (2) 

where 4>' is called the dual space of $ and contains the linear functionals over 4>. 

The basic reason why we need the spaces 4>' and 4> x is that the bras and kets associ- 
ated with the elements in the continuous spectrum of an observable belong, respectively, 
to <J>' and 3» x rather than to M '. The basic reason reason why we need the space 4> is that 
unbounded operators are not defined on the whole of but only on dense subdomains 
of that are not invariant under the action of the observables. Such non-invariance 
makes expectation values, uncertainties and commutation relations not well defined on 
the whole of Jrff. The space 4> is the largest subspace of the Hilbert space on which such 
expectation values, uncertainties and commutation relations are well defined. 

Besides accommodating resonances and Dirac's bra-ket formalism, the rigged Hilbert 
space seems to capture the physical principles of quantum mechanics better than the 
Hilbert space. For example, assuming that the Hilbert space provides the whole math- 
ematical framework for quantum mechanics leads to the conclusion that Heisenberg's 
uncertainty relations are not physical, since they cannot be defined on the whole of the 
Hilbert space [1]. Using the rigged Hilbert space, one overcomes this difficulty after 
realizing that the commutation relations are well defined on 4>. 

The completeness relation is a good place to appreciate the added value of the rigged 
Hilbert space. Consider, for example, the Hamiltonian H of a system. In the Hilbert 
space, one writes the completeness relation as 

1= / dE E , (3) 

where Eg are the spectral projections of H and Sp(H) is its spectrum. However, within 
the rigged Hilbert space one can write 1 

poo 

1 = £|E M )(E M |+ / dE\E)(E\, (4) 

where \E n ) and \E) are the bound and scattering states of H, respectively. In addition 
to (4), the rigged Hilbert space gives you an additional completeness relation in which 
the resonance states participate: 

l=£|E M )(E M |+£|z n )(z n |+ / dE\E){E\, (5) 

where \z n ) are the Gamow (resonance) states of H and the last integral, called the back- 
ground, is performed in the complex plane right below the negative real axis of the 



In the Hilbert space, one can actually write something close to, although not the same as (4), by means 
of direct integral decompositions. 
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second sheet. 2 Thus, the completeness relation (5) substitutes the scattering states con- 
tribution by the resonance contribution plus a background, thereby putting the resonance 
spectrum on the same footing as the bound and scattering spectra. 

It is important to note that the integrals in (4) and (5) are different, and that the 
resonance contribution does not appear in (4), because resonances are not asymptotic 
states. Also important is to note that the resonance states, and therefore expansion (5), 
need a different rigged Hilbert space from that needed by the scattering states and 
expansion (4). 

There are dangers in using the rigged Hilbert space indiscriminately, though. For in- 
stance, A. Bohm and collaborators have been using a rigged Hilbert space (of Hardy 
class) to construct a quantum theory of resonances, see review [2] and references 
therein. However, such theory is inconsistent with quantum mechanics and must be dis- 
carded [3]. 

The structure of these lectures is as follows. After an introductory Lecture 1, I will 
explain in Lecture 2 how to construct the rigged Hilbert space of the one-dimensional 
rectangular barrier potential. This will show that the rigged Hilbert space is already 
needed to provide the mathematical support of the most basic quantum systems. In 
Lecture 3, 1 will construct the rigged Hilbert space of the Lippmann-Schwinger equation, 
which equation governs quantum scattering. In Lecture 4, I will construct the rigged 
Hilbert space of the analytic continuation of the Lippmann-Schwinger equation. The 
rigged Hilbert space of Lecture 4 will be needed in the final Lecture 5 to provide the 
mathematical support for the resonance states. The PDF file of each talk will be posted 
at http : / / www .ucsd.edu/~rafa. 

Since space prevents a full account, at each lecture I will refer the reader to the 
appropriate papers where further details can be found. The essentials of functional 
analysis needed to understand those papers can be found in [4, Chapter 2]. 

CONSTRUCTION OF A SIMPLE RIGGED HILBERT SPACE: 

LECTURE 2 

If resonances are not included, the way to construct the rigged Hilbert space of a 
quantum system is as follows: 

1. We identify the observables of the system. Their expressions are usually given by 
linear, differential operators. 

2. We identify the Hilbert space, whose scalar product is used to calculate probability 
amplitudes. 

3. We identify the domains, spectra and eigenfunctions of the observables. If the 
observables have a discrete spectrum, we need not go beyond the Hilbert space. 
However, if at least one of the observables has continuous spectrum, we need to 
enlarge the Hilbert space to the rigged Hilbert space. (If position and momentum 



2 As we will see in Lecture 5, the expansion (5) must be either regulated or understood in a time-dependent 
way. 
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where 



are among the observables, we will always need the rigged Hilbert space.) 

4. We construct the space 4> in which physical quantities such as expectation values, 
uncertainties and commutation relations are well defined. When resonances are not 
included, the space 4> is usually given by the maximal invariant subspace of the 
algebra of observables. 

5. We construct the dual <J>' and antidual 4> x spaces. We construct the bras and kets of 
the observables and check that they respectively belong to 4>' and 3» x . 

6. The completeness relations and all the features of Dirac's bra-ket formalism now 
follow. 

Let's see how the above steps are carried out in the case of a spinless particle moving 
in one dimension and impinging on a rectangular barrier. The observables relevant to 
this system are the position Q, the momentum P, and the Hamiltonian H: 

Qf(x)=xf(x), (6) 
Pf( x ) = -ih-jLf(x), (7) 

— oo < x < a 

V(x) = { Vq a<x<b (9) 
b < x < oo 

is the ID rectangular barrier potential. These observables satisfy the following commu- 
tation relations: 

[Q,P] = ihI, (10) 
ih 

[H,Q] = P, (11) 

m 

dV 

[H,P] = ih—. (12) 

Since our particle can move in the full real line, the Hilbert space on which the 
differential operators (6)-(8) should act is 

/oo 
dx\f(x)\ 2 <oo}. (13) 

The corresponding scalar product is 

/CO 
dxf(x)g(x), f,geL 2 . (14) 
-oo 

The differential operators (6)-(8) induce three linear operators on the Hilbert space 
L 2 . These operators cannot be defined on the whole of L 2 , but only on the following 
subdomains of L 2 : 

@(Q) = {feL 2 \xfeL 2 }, (15) 
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®{P) = {feL 2 \feAC,Pfe L 2 } , (16) 
0(H) = {f e L 2 \ f e AC 2 , Hfel 2 }, (17) 

where, essentially, AC is the space of functions whose derivative exists, and AC 2 is the 
space of functions whose second derivative exists. On these domains, the operators Q, P 
and H are self-adjoint, and their spectra are 

Sp(<2) = Sp(P) = (-00,00), Sp(tf) = [0,oo), (18) 

which spectra coincide with those we would expect on physical grounds. 

To obtain the eigenfunctions corresponding to each eigenvalue, we have to solve the 
eigenvalue equation for each observable: 

x(x\x') = x' (x\x') , (19) 
-i^(x\p) = p(x\p) , (20) 

("l;i +F w) WE>=£W£> - <21) 

The eigenfunctions of Q are delta functions, 

(x\x / ) = 8(x-x), (22) 



those of P are plane waves, 



and those of H are given by 



Jpx/h 



, m xi/2 f T(k)e- lk * -™<x<a 

(x\E + ) T =l 2 xl Ar{k) iKX + B t (k)e- iKX a<x<b (24) 

V 2nkh J y R r (k)e ikx + e- ikx b<x<oo : 



, xl/2 f e lkx + Ri(k)e- lkx -<*><x<a 
(x\E + ) x = ( — J x <^ Ai(*)e iKX + Bi(ik)e- iKX a < x < b 
\2nkh J I T{k)e ikx b<x<oo, 



(25) 



where 



k= X / TT E, k=J- t (E-V ), (26) 



' 2m / 2m 

and where the coefficients that appear in Eqs. (24)-(25) can be easily found by the 
standard matching conditions at the discontinuities of the potential. Physically, (x\E + ) T 
((x\E + )i) represents a particle of energy E impinging on the barrier from the right (left). 

The eigenfunctions (22)-(25) are not square integrable, that is, they do not belong to 
L 2 . Mathematically speaking, this is the reason why they are to be dealt with as distribu- 
tions (note that all of them except for the delta function are also proper functions). 
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We now start the construction of the rigged Hilbert space by constructing 4>. The 
space 4> is given by 

oo 

cj>= p| 3>(A n B m ). (27) 

n,m—0 
A,B=Q,P,H 

In view of expressions (6)-(8), 4> is simply 

$ = {(pGL 2 | <peC°°(M), <p( n \a) = q>W(b) = 0, n = 0,l,..., 

P n Q m H l (p(x) GL 2 , n,m,/ = 0,l,...}, (28) 

where C°°(R) is the collection of infinitely differentiable functions, and <pM denotes the 
nth derivative of <p. From the last condition in Eq. (28), we deduce that the elements of 
4> satisfy the following estimates: 



\\(p\\n,m,l = \j J_ dx\P n Q m W(p{x)\ 2 <oo, „,!»,/ = 0, 1, ... . (29) 

These estimates mean that the action of any combination of any power of the observables 
remains square integrable. For this to happen, the functions <p(x) must be infinitely 
differentiable and must fall off at infinity faster than any polynomial. Hence, 4> is a 
Schwartz-like space. 

Because 4> is invariant under the action of the observables, 



A<fC<D, A = P,Q,H, (30) 

the expectation values 

(<M» , <P e <*>, A = P, Q,H, n = 0, 1, . . . (31) 

are finite, and the commutation relations (10)-(12) are well defined. In particular, Heisen- 
berg's uncertainty principle makes sense on 4>. 

The spaces <J>' and 4> x are simply the collection of linear and antilinear functionals 
over 4>, respectively. By combining the spaces <f>, J^f, 4> x and <J>', we obtain the rigged 
Hilbert spaces of our system, 

<J>C^rc4> x , (32) 

<DcJc«J>'. (33) 

The space 3> x accommodates the eigenkets \p), \x) and \E + )\ S of P, Q and H, whereas 
4>' accommodates the eigenbras (p\, (x\ and \ r { + E\. 

Mathematically, the bras and kets are distributions defined as follows. Given a func- 
tion f(x) and a space of test functions the antilinear functional F that corresponds to 
the function f(x) is an integral operator whose kernel is precisely f(x) : 

F((p) = Jdxtfx)f(x), (34) 



6 



and the linear functional F generated by the function f(x) is an integral operator whose 
kernel is the complex conjugate of f(x): 

F{<p) = Jdx(p(x)f(x). (35) 

In Dirac's notation, these two equations become 

(<p\F) = J dx(<p\x)(x\f), (36) 

(F\<p) = fdx(f\x)(x\<p). (37) 

Note that these definitions are very similar, except that the complex conjugation affects 
either f(x) or (p(x), which makes the corresponding functional either linear or antilinear. 

Definitions (34) and (35) provide the link between the quantum mechanical formalism 
and the theory of distributions. In practical applications, what one obtains from the 
quantum mechanical formalism is the distribution f(x) (in this lecture, the plane waves 
-j=e ipx / h , the delta function 8(x — x') and the eigenfunctions (x\E + )i jT ). Once f(x) is 

given, one can use definitions (34) and (35) to generate the functionals |F) and (F\. 
Then, the theory of distributions can be used to obtain the properties of the functionals 
\F) and (F\, which in turn yield the properties of the distribution f(x). 

By using prescription (34), we can define for each eigenvalue p the eigenket \p) 
associated with the eigenfunction (23): 



(q>\p) = j_Jxcp(x)-j=e^/\ (38) 
which, using Dirac's notation for the integrand, becomes 

/oo 
dx(<p\x)(x\p). (39) 
-oo 

Similarly, for each x, we can define the ket |jc) as 

/oo 
dx (p(x')8(x-x), (40) 
-co 

which, using Dirac's notation for the integrand, becomes 

/oo 
dx ((p\x)(x\x). (41) 
-co 

The definition of the kets |2s + )i r that correspond to the Hamiltonian's eigenfunc- 
tions (24)-(25) follows the same prescription: 

/oo 
dx^(x)(x\E + \ r , (42) 
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that is, 

/oo 
dx(<p\x)(x\E+) hr . (43) 

One can now show that the definition of the kets \p), \x) and \E + )\ S makes sense, and 
that these kets indeed belong to the space of distributions 4> x . 

By using prescription (35), we can also define for each eigenvalue p the eigenbra (p\ 
associated with the eigenfunction (23): 

/oo 1 poo 

Jx 9 {x)-^r^l h = j_Jx(p\x)(x\(p). (44) 

Comparison with Eq. (38) shows that the action of (p\ is the complex conjugate of the 
action of \p), 

{p\<p)=Jv\pj, (45) 

and that 

(p\x) = Wfi = -±=e-W h . (46) 

The bra (x\ is defined as 

/oo poo 
dx (p(x)8(x-x) = / dx (x\x)(x\(p). (47) 
-CO J — CO 

Comparison with Eq. (40) shows that the action of (x| is complex conjugated to the 
action of \x), 

(x\(p)=Wj, (48) 

and that 

(x\x) = {x\x} = 8(x-x). (49) 
Analogously, the eigenbras of the Hamiltonian are defined as 

/oo p CO 

dx(p(x) hr ( + E\x)= dx l!T ( + E\x)(x\(p), (50) 
-oo J — oo 

where 



Kr ( + E\x) = (x\E+) lr (51) 

Comparison of Eq. (50) with Eq. (42) shows that the actions of the bras i, r ( + £| are the 
complex conjugates of the actions of the kets \E + )iy. 



Kr ( + E\( P ) = (( P \E+) ljI . (52) 

Now, by using the rigged Hilbert space, one can show that the definitions of (p\, (x\ and 
l,r{ + E\ make sense and that (p\, (x\ and i, r ( + £'| belong to 4>'. 

It is important to keep in mind the difference between eigenfunctions and kets. For 
instance, (x\p) is an eigenfunction of a differential equation, Eq. (20), whereas \p) is a 
functional, the relation between them being given by Eq. (39). A similar relation holds 



between (x'\x) and \x), and between (x\E + )\ s and \E + )\j. It is also important to keep in 
mind that "scalar products" like (x\p), (x!\x) or {x\E + )^ r do not represent an actual scalar 
product of two functionals; these "scalar products" are simply solutions to differential 
equations. 

The kets \p), \x) and \E + )i jT are indeed eigenvectors of P, Q and H, respectively: 

P\p)=p\p), peR, (53) 

Q|jc) = jc|jc) , xGM, (54) 
H\E+) hl = E\E + ) l ^ £€[0,00). (55) 
Similarly, the bras (p\, (x\ and i ;r ( + £'| are left eigenvectors of P, Q and H, respectively: 

(p\P = p{p\, (56) 

(x\Q = x(x\, xGK, (57) 

lyT { + E\H = E 1 , I {+E\, £e[0,oo). (58) 

Note that these equations are to be understood in the distributional way, that is, as 
"sandwiches" with elements of 4>. For example, Eq. (53) should be understood as 

(<p\P\p) = p(<P\p) , peR^ye®, (59) 

and the same for (54)-(58). Usually, the "sandwiching" is implicit and therefore omitted. 

Now that we have constructed the Dirac bras and kets, we can see how other aspects 
of Dirac's bra-ket formalism hold within the rigged Hilbert space. For example, the 
completeness relations 

poo 

dp\p)(p\=I, (60) 



f 



dx'\x')(x'\ =/, (61) 



/'OO /' oo 

/ dE\E + } u {+E\+ dE\E + ) II {+E\=I, (62) 
Jo JO 

and the action of P, Q and H, 



P= / dpp\p)(p\, (63) 

J co 

/CO 
dxx\x) (x\ , (64) 
-co 



H= / dEE\E + ) u ( + E\ + / dE E\E + ) TI ( + E\ , (65) 
Jo Jo 

all hold within the rigged Hilbert space as a "sandwich" with elements of Other 
expressions such as the delta normalization of eigenfunctions, e.g., 

1 /'CO 

^J_J™ i{p - p ' )x/h = S(p-p')i (66) 
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or the "matrix elements" of the observables, e.g., 

(x\Q\x) =x 8(x-x), (67) 

(x\P\x f ) = -ih4- 8(x-x r ), (68) 
ax 

M"M = (-|;^m)^-A (69) 

are interpreted the same way. 

To conclude this section, I would like to refer the reader to [5, 6] for a detailed account 
of the above results. 



THE RIGGED HILBERT SPACE OF THE 
LIPPMANN-SCHWINGER EQUATION: LECTURE 3 

The Lippmann-Schwinger equation is one of the cornerstones of scattering theory. It is 
written as 

^ E±) ^ E)+ E^k^ nE±) ' <70) 

where jE*) are the "in" and "out" Lippmann-Schwinger kets, \E) is an eigenket of the 
free Hamiltonian Hq, 

H \E)=E\E), (71) 

and V is the potential. The Lippmann-Schwinger kets are, in particular, eigenvectors of 
//: 

H\E ± )=E\E ± ). (72) 
To the kets |£' ± ), there correspond the bras ( ± £'|, which satisfy 

(±£| = < £|+<±£|v £d^ <73) 

The bras ( ± £'| are left eigenvectors of H, 

{ ± E\H = E{ ± E\, (74) 

and the bras (E\ are left eigenvectors of Hq, 

(E\H = E(E\. (75) 

The Lippmann-Schwinger equation (70) for the "in" \E + ) and "out" kets has 
the scattering "in" and "out" boundary conditions built into the ±ie, since Eq. (70) 
is equivalent to the time-independent Schrodinger equation (72) subject to those "in" 
(+ie) and "out" (— ie) boundary conditions. In the position representation, the ±i£ 
prescriptions yield the following asymptotic behaviors: 



(x\E + )-~ >e lkz + f(k,G) e -^, (76) 
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(x\E-) ► e lkz + f(k, 9) , (77) 

y — foo f 

where x = (x,y,z) = (r, 0,0) are the position coordinates, k is the wave number of 
Eq. (26) and f(k, 9) is the scattering amplitude. 

Likewise any bra and ket, the Lippmann-Schwinger bras and kets do not have a place 
in the Hilbert space. In this lecture, we will construct the rigged Hilbert space to which 
they belong. We will use the example of the spherical shell potential. 



The radial Lippmann-Schwinger equation 

For the spherical shell potential, 

{0 0<r<a 
V a<r<b (78) 
b<r 

the Lippmann-Schwinger equation can be solved explicitly. Because the potential (78) 
is spherically symmetric, we will work in the radial position representation and restrict 
ourselves to angular momentum / = 0. 

In the radial representation, the Lippmann-Schwinger equation (70) becomes 

{r\E±) = (r\E) + ( r | — J-— - V|£±) . (79) 

The procedure to solve Eq. (79) is well known. Since Eq. (79) is an integral equation, it is 
equivalent to a differential equation subject to the boundary conditions that are built into 
it. In our case, for / = 0, Eq. (79) is equivalent to the Schrodinger differential equation, 

■££ + V(r)) W *±)=E<r|E±>, (80) 

subject to the following boundary conditions: 

{r\E ± ) = 0, (81) 
(r|£ ,± ) is continuous at r = a, b , (82) 

-^-(r|£ ,± ) is continuous at r = a, b, (83) 
dr 

(r\E+)~ e - ikr -S(E)e ikr asr^oo, (84) 

(r\E~) ~e ikr -S(E)e- ikr asr^oo, (85) 

where S(E) is the S matrix in the energy representation. The boundary conditions (84) 
and (85) originate from the ±/e conditions of Eq. (79). The asymptotic behaviors (84) 
and (85) are the / = 0, radial counterparts of the asymptotic behaviors (76) and (77). 
If we insert (78) into (80), and solve (80) subject to (81)-(85), we obtain 

(rl^) ee xHnE) = N(E) ^f^- , E e [0, °°) , (86) 
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where N(E) is a delta-normalization factor, 

N{E) = 

\\ n sj2m/h 2 E 
X(f",E) is the so-called regular solution of Eq. (80), 



1 2m/K (87) 



sin( A /|££r) 0<r<a 

2m ( c \7 \ : 2m 



X{r,E)=l M E)e i ^ {E - Vo)r +ME)e- i ^ {E - V ° )r a<r<b TO 



J z {E)e^ r + J A {E)e W"' b<r<°°. 



and ^±{E) are the Jost functions, 

J?+(E) = -2i t / 4 (E), (89) 

e /-(E)=2i e / 3 (E). (90) 

The explicit expressions for J?<x can be obtained by matching the values of %{r;E) 
and of its derivative at the discontinuities of the potential. In terms of the Jost functions, 
the S matrix is given by 

Similarly to the "right" ones, we can obtain the "left" Lippmann-Schwinger eigen- 
functions by writing Eq. (73) in the radial position representation, 

(±E | r ) = (£ | r ) + (±E | Fl -±_-j r ). (92) 

Solving Eq. (92) is analogous to solving Eq. (79). The solutions to (92) are the complex 
conjugates of the solutions to (79): 



( ± E\r)= X ± (nE)=X T (r;E). (93) 



The rigged Hilbert space of the Lippmann-Schwinger bras and kets 

By using (34), we associate "in" and "out" kets |£' ± ) with the eigenfunctions x ± { r ',E) 
for each E e [0, °°) : 

POO POO 

(^1^)=/ dryHr)x ± (r,E)= dr (<p±|r> (r^) , (94) 
Jo Jo 

where (p^ belong to the space 4> that will be constructed below. (Note that the 4> of this 
lecture is different from the 4> of Lecture 2.) Similarly, by using (35), we define the bras 

{ ± E\ as 

poo POO 

(±£1^)=/ drx ± (r,E)(p ± (r)= / drPEMMty*) . (95) 
Jo Jo 
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Note that even though % (r;E) = (r\E ) are also meaningful for complex energies, the 
energy in Eqs. (94) and (95) runs only over Sp(H) = [0,°°), because in this lecture we 
restrict ourselves to bras and kets associated with energies that belong to the scattering 
spectrum of the Hamiltonian. 

From definitions (94) and (95), it follows that the action of ( ± £'| is complex conju- 
gated to the action of {E^: 

( ± £|9 ± ) = (9±]£±). (96) 

We now need to find the subspace 4> on which the above definitions make sense. 
Besides making (94)-(95) well defined, the space 4> must also be invariant under the 
action of the observables of the system. Since in this lecture the only observable we are 
concerned with is the Hamiltonian, we will simply require invariance under H. Thus, the 
space $ must satisfy the following conditions: 

• The space <J> is invariant under the action of H. (97) 

• The elements of 4> are such that the integrals in Eqs. (94) -(95) make sense. (98) 

In order to meet requirement (97), the wave functions <p ± (r) must at least be in the 
maximal invariant subspace of H: 

CO 

& = p| Q){H n ) . (99) 

n=0 

In order to meet requirement (98), the wave functions <p ± (r) must behave well enough 
so the integrals in Eqs. (94)-(95) are well defined. From the expression for % ± (r;E), 
Eq. (86), one can see that the ^) ± (r) have essentially to control purely imaginary expo- 
nentials. Therefore, the space $ that meets the requirements (97)-(98) is 

cD={(p± G L 2 ([0,oo),Jr)| <p±e^, ||(p ± ||„ ;ra <oo, n,m = 0,1,2, ■■■}, (100) 

where the || || M;OT are given by 



||9 ± ll»^ = Y < / ^|(l+r)"(l+//)™<p±(r)| 2 , n,m = 0,l,2,.... (101) 

The space 4> is thus the collection of square integrable functions that belong to the max- 
imal invariant subspace of H and for which the estimates (101) are finite. In particular, 
because <p ± (r) satisfy the estimates (101), <p ± (r) fall off at infinity faster than any poly- 
nomial of r. 

lim(l + r) n <p ± (r) = 0, n = 0,1,2,.... (102) 

r — >oo 

Thus, likewise in Lecture 2, we have obtained a Schwartz-like space. Obviously, the 
space $ can also be seen as the maximal invariant subspace of the algebra generated by 
the Hamiltonian and the operator multiplication by r. 

Note that we have used superscripts ± to denote the elements of one and the same 
space The reason is that when we use +, it means that the elements of 4> are acted 
upon by \E + ) or { + E\, and when we use — , they are acted upon by or 
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We can now construct the spaces and 4> x , and see that lE^) belong to <3> x , whereas 
( ± £'| belong to <&'. Thus, the Lippmann-Schwinger equations (70) and (73), Eqs. (72) 
and (74), as well as 

e -mt/h\ E ±} = g -iEt/h\ E ±} ^ (103) 

(± E ^ e -iHt/h = e iEt/h ( ± E ^ (104) 

hold in the distributional sense, i.e., as "sandwiches" with elements of 4>. 

For a detailed account on the results of this lecture, the reader may wish to refer to [7]. 

THE RIGGED HILBERT SPACE OF THE ANALYTIC 
CONTINUATION OF THE THE LIPPMANN-SCHWINGER 
EQUATION: LECTURE 4 

This lecture is devoted to construct and characterize the analytic continuation of the 
Lippmann-Schwinger bras and kets. As in Lecture 3, we restrict ourselves to the spheri- 
cal shell potential (78) and zero angular momentum. 

The ultimate goal we want to achieve by analytically continuing the solutions of 
the Lippmann-Schwinger equation is to obtain, in Lecture 5, the resonance (decay) 
amplitude. 



The wave number representation 

The Lippmann-Schwinger eigenfunctions depend explicitly on the wave number k 
of Eq. (26) rather than on the energy E. It is therefore convenient to rewrite their 
expressions in terms of k before performing analytic continuations. 

We start by writing the regular solution (88) in terms of k: 

{sin(fcr) < r < a 

Mk)e iKr + Mfy~ iKr a<r<b (105) 
J? 3 (k) e ikr + f A {k)e- ikr b < r < °o , 

where K is given by (26). In terms of k, the Lippmann-Schwinger eigenfunctions read as 



The eigenfunctions % ± (r;E) are 5-normalized as functions of E. The Lippmann- 
Schwinger eigenfunctions that are 5-normalized as functions of k are given by 



fM) S ^^ M =^g. (107) 

In bra-ket notation, we will write 

(r\k ± } = x ± (r,k), k>0, (108) 
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( ± k\r) = x ± (nk)=x T (nk) 1 k>0. (109) 
Because of (107), in terms of k the Lippmann-Schwinger bras and kets read as 



( ± k\ = \l ^-2k( ± E\, k>0, (110) 
" 2m 



\k ± ) = \\^-2k\E ± ), k>0. (Ill) 

The bras ( ± fc| and kets \k^) are, respectively, left and right eigenvectors of H with 
eigenvalue j^k 2 : 

h 2 

(±k\H=—k 2 ( ± k\, (112) 



2m 
h 2 

H\k ± ) = — k 2 ^). (113) 
2m 



The analytic continuation of the Lippmann-Schwinger eigenfunctions 

The analytic continuation of X ± ( r >k) is done as follows. First, one specifies the 
boundary values that the Lippmann-Schwinger eigenfunctions take on the positive k- 
axis. And second, one continues those boundary values into the whole fc-plane. Since 
the boundary values of the Lippmann-Schwinger eigenfunctions on the positive &-axis 
are given by Eq. (107), and since x ± (f",k) are expressed in terms of well-known analytic 
functions, the continuation of % ± (r;k) from the positive fc-axis into the whole wave- 
number plane is well defined. 

A word on notation. Whenever they become complex, we will denote the energy 
E and the wave number k by respectively z and q. Accordingly, the continuations of 
X ± (r;E) and x ± (f',k) will be denoted by X ± ( r >z) and X ± ( r >q)- m bra-ket notation, the 
analytically continued eigenfunctions will be written as 

(r\q ± )=X ± (r,q), (114) 
( ± q\r) = X T (r,q). (115) 

Note that, in distinction to (93), the analytically continued "left" eigenfunction is not the 
complex conjugate of the "right" eigenfunction but 

( ± q\r) = {r\q t ). (116) 
In order to characterize the analytic properties of ^ ± (r;^), it is useful to define 

Z ± = {qeC\ jr±(q)=0}. (117) 

The elements of Z + are simply the resonance energies. Since xi^q) an d c /±(q) are 
analytic in the whole fc-plane, X ± ( r ->q) is analytic in the whole fc-plane except at Z±, 
where its poles are located. 
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In order to define the analytically continued Lippmann-Schwinger bras and kets, we 
need to know how X ± ( r > < l) grow with q. To find out, let us recall first that the growth of 
%{r\q) is bounded by [8] 

|z(r^)|<C-^-el Im ^l^, qeC. (118) 

From Eqs. (107) and (118), it follows that the eigenfunctions % ± (r;q) satisfy 

\x^)\<Y^T¥fr ellm{q)K < 119) 
1 1 \f±{q)\ i + kk 

When q e Z±, % ± (r;q) blows up to infinity. 

We can further refine the estimates (1 19) by characterizing the growth of l/\ c f±(q) | 
in different regions of the complex plane. In the upper half of the complex fc-plane, the 
inverse of J?+(q) is bounded: 

<C, \m(q)>0. (120) 



I ,/+(*)! 

l 



In the lower half -plane, is infinite whenever qeZ + . As \q\ tends to °° in the lower 

half plane, we have 

\q\^~,lm(q)<0). (121) 



f+{q) \-Cq- 2 e 2i « bl 

The above estimates are satisfied by (q) when we exchange the upper for the lower 
half plane, and Z + for Z . 



The analytic continuation of the Lippmann-Schwinger bras and kets 

The analytic continuation of the Lippmann-Schwinger bras is defined for any complex 
wave number q in the distributional way (35): 

poo noo 

(^1^)=/ dr(p ± (r) X T (r,q) = dr (± q\r) (r^) , (122) 
JO JO 

where the functions <p ± (r) belong to a space of test functions 4> exp that will be 
constructed below. Similarly to the bras, the analytic continuation of the Lippmann- 
Schwinger kets is defined by way of (34): 

<p=V>=/ dr~^(rjx ± (r,q)= dr {<p±\r) {r\q±) . (123) 
JO JO 

Note that definition (122) is actually a slight generalization of (35), due to (1 16). 

The bras (122) and kets (123) are defined for all complex q except at those q at which 
the corresponding eigenfunction has a pole. At those poles, one can still define bras and 
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kets if in definitions (122) and (123) one substitutes the eigenfunctions X ± ( r > ( l) by their 
residues at the pole. 

From the analytic continuation of the bras and kets into any complex wave number, 
one can now obtain the analytic continuation of the bras and kets into any complex 
energy of the Riemann surface (compare with Eqs. (1 10) and (111)): 




1 2m 1 



\z ± ) = \l — -\<l ± )i ( ±z \ = \l-^Y q {±ql (124) 



Construction of the rigged Hilbert space 

Likewise the bras and kets associated with real energies, the analytic continuation 
of the Lippmann-Schwinger bras and kets must be described within the rigged Hilbert 
space rather than just within the Hilbert space. We will denote the rigged Hilbert space 
for the analytically continued bras by 

^>expCL 2 ([0,oo),Jr)c^ xp , (125) 

and the one for the analytically continued kets by 

*ex P CL 2 ([0,oo),Jr)c4> e x xp . (126) 

The functions ^ e 4> exp must satisfy the following conditions: 

• They belong to the maximal invariant subspace Q) of H, see Eq. (99). (127) 



They are such that definitions (122) and (123) make sense. 



(128) 



The reason why <p ± must satisfy condition (127) is that such condition guarantees that 
all the powers of the Hamiltonian are well defined. Condition (127), however, is not 
sufficient to obtain well-defined bras and kets associated with complex wave numbers. 
In order for ( ± ^r| and \q^) to be well defined, the wave functions <p ± (r) must be well 
behaved so the integrals in Eqs. (122) and (123) converge. Since by Eq. (119) X ± ( r >4) 
grow exponentially with r, the wave functions ?> ± (r) have to, essentially, tame real 
exponentials. If we define 



dr 



nr 



l+nr 



n,n = 0, 1,2,. 



(129) 



then the space 4> exp is given by 



exp 



{(p ± e$\ ||fl> ± ||„ )B /<oo, «y = 0,1,2,...} 



(130) 



This is just the space of square integrable functions which belong to the maximal 
invariant subspace of H and for which the quantities (129) are finite. In particular, 
because <p ± (r) satisfy the estimates (129), their tails fall off faster than Gaussians. 



17 



From Eq. (119), it is clear that the integrals in Eqs. (122) and (123) converge already 
for functions that fall off at infinity faster than any exponential. We have imposed 
Gaussian falloff because it will allow us to perform resonance expansions in Lecture 5. 

It is illuminating to compare the space <f> of Lecture 3 with the space 4> exp of Eq. ( 1 30). 
Because for real wave numbers the Lippmann-Schwinger eigenfunctions behave like 
purely imaginary exponentials, in Lecture 3 we only needed to impose on the test 
functions a polynomial falloff, thereby obtaining a space of test functions very similar to 
the Schwartz space. By contrast, for complex wave numbers the Lippmann-Schwinger 
eigenfunctions blow up exponentially, and therefore we need to impose on the test 
functions an exponential falloff that damps such an exponential blowup. 

One can now easily show that the kets {q^ belong to <I>* and satisfy 



H\q ± ) = ^-q 2 \q ± ), (131) 

h 2 

( ± q\H=—q 2 ( ± q\, (133) 

^qle^'^ = e ^ 2fe /( 2m ) (± ? | . (134) 
Equations (131) and (133) can be rewritten in terms of the complex energy z as 

H\z ± )=z\z ± ), (135) 
{ ± z\H = z{ ± z\. (136) 

Note that (136) is not given by ( ± z|// = z( ± z|, as one may naively expect from formally 
obtaining (136) by Hermitian conjugation of (135). 

For a full account of this lecture, the reader can refer to [9]. 



Similarly, ( ± q\ e <J>g Xp and 



RESONANCE STATES, AND THEIR RIGGED HILBERT SPACE: 

LECTURE 5 

The Gamow states are the state vectors of resonances. They are eigenvectors of the 
Hamiltonian with a complex eigenvalue. The real (imaginary) part of the complex 
eigenvalue is associated with the energy (width) of the resonance. 

Because self-adjoint operators on a Hilbert space can only have real eigenvalues, the 
Gamow states fit not within the Hilbert space but within the rigged Hilbert space. In 
this lecture, we will see that the Gamow states belong to the rigged Hilbert space of 
Lecture 4. 

Like in Lectures 3 and 4, we will use the spherical shell potential (78) and study the 
zero angular momentum case only. Unlike in Lecture 4, we will write most results in 
terms of the energy, because they tend to be simpler than in terms of the wave number. 
The energy and the wave number of a resonance R will be denoted by zr and &r. 
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The Gamow eigenfunctions satisfy the Schrodinger equation 

h 2 d 2 



Imdr 2 +V ^J M ( r;ZR ) =ZRu(r;zR), (137) 
subject to "purely outgoing boundary conditions," 

w(0;z R )=0, (138) 

u(r;zR) is continuous at r = a,b , (139) 
d 

—u(r;zR) is continuous at r = a,b , (140) 
dr 

u(r; Z R) ~e* Rr asr^oo, (141) 

where (141) is the "purely outgoing boundary condition" (POBC). Comparison of (138)- 
(141) with (81)-(85) shows that it is the POBC what selects the resonance energies. 

For the potential (78), the only possible eigenvalues of (137) subject to (138)-(141) 
are the zeros of the Jost function, 

< /+(zr)=0. (142) 

The solutions of this equation come as a denumerable number of complex conjugate 
pairs Zni The number Zn = E n — iT„/2 is the nth resonance energy, and z* = E n + iT n /2 
is the nth anti-resonance energy. The corresponding wave numbers are 



/ 2m / 2m 

k n = d-^Zn, -k* n = J T2~4; /i =1,2,.... (143) 

For the potential (78), the resonance energies are simple poles of the S matrix (see [10] 
for an example of a potential that produces double poles). In order to write expressions 
for resonances and anti-resonances together, we will label the resonances by a positive 

integer n = 1,2,... and the anti-resonances by a negative integer n = — 1 , —2, 

The nth Gamow eigensolution, n = ±1, ±2, . . ., reads 



73M sin (M 0<r<a 

Mkn) 

? ik n r b < r < 



u{r;zn) = u{r;k n )=N n { j^e iQ " r + ^ j e~ iQnr a<r<b (144) 



where 



/ 2m 

Qn=\/-^(Zn-V ), (145) 



and N n is a normalization factor, 

t2 



N l n =i^[S(q)] q=kn . (146) 
The Gamow eigenfunctions (144) are related with % ± (r;q) by 



v 2tt 

u(r,k n ) = -^ res [X + {r;q)] q=kn , (147) 
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u(r;k n ) = iV27uN n X~(r,k n ) . (148) 
Because of (148) and (116), the "left" Gamow eigenfunctions read 

(z n \r) = [u(r;z* n )}*, * = ±1,±2,. .. . (149) 

Thus, the "left" Gamow eigenfunction is not just the complex conjugate of the "right" 
eigenf unction, but the complex conjugated eigenfunction evaluated at the complex con- 
jugated energy. Because the Gamow eigenfunctions satisfy 

[u(r;z* n )}* = u(r;z n ), n = ±l,±2,..., (150) 

the "left" and the "right" Gamow eigenfunctions are actually the same eigenfunction, 

(z n \r) = [u(r;z* n )}* = u(r;z n ) = (r\z n ) , n = ±1,±2, . . . . (151) 

In terms of the wave number, Eq. (151) reads as 

(kn\r) = [u(r--kl)Y =u(r;k n ) = (r\k n ) , n = ±1,±2,... . (152) 

The Gamow eigenfunctions u(r;z n ) are not square integrable and therefore must be 
treated as distributions. By treating them as distributions, we will be able to generate 
the Gamow bras and kets. According to (34), the Gamow ket \z n ) associated with the 
eigenfunction u(r;z n ) must be defined as 

POO poo 

{<p\z n )= dr[<p(r)}*u(r,z n )= dr (<p\r) (r\z n ) , n = ±l,±2,.... (153) 
Jo Jo 

Similarly, the Gamow bra associated with the resonance (or anti-resonance) energy z n is 
defined as 

POO POO 

(z n \(p)= dr(p(r)u(r;z n )= dr (z n \r) (r\(p) , n = ±l,±2,.... (154) 
JO JO 

From these definitions and from Eqs. (147) and (148), it is clear that the Gamow bras 
and kets are accommodated by the rigged Hilbert spaces (125) and (126). 

Within the rigged Hilbert spaces (125) and (126), it holds that the Gamow bras and 
kets are eigenvectors of the Hamiltonian: 

H\z n ) = Z n \z n ) , »= ±1,±2, .. . , (155) 

(z n \H = z n {z n \ , n = ±l,±2,.... (156) 

Because their energy is complex, the time evolution of the Gamow bras and kets should 
be time asymmetric. For resonances, it should be that 

(z n \e- iHt / h = e iZnt / h (z n \, onlyfor?<0, n=l,2,..., (157) 

e - iHt / h \z n ) = e - iZnt l n \z n ), onlyfor?>0, n=l,2,..., (158) 
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whereas for anti-resonances, it should be that 



(z n \e- iHt/h = e iZnt l h {z n \ , only for t > 0, n = -1, -2, . . . , (159) 

e - iHt/h \z n ) = e- iZnt/h \z n ) , only for t < 0, n = -1, -2, . . . . (160) 
If we define the complex delta function at z n by 

dEf(E)S(E-Zn)=f(Zn), (161) 



Jo 



one can use the results of Lecture 4 to show that the Gamow eigenfunction u(r;z n ), 
the complex delta function (multiplied by a normalization factor) and the Breit-Wigner 
amplitude (multiplied by a normalization factor) are linked with each other: 

u{r;z n ) <-> iV27cJr„8(E-z n ), Ee [0,oo) <-> , £ e (-00,00) 

posit, repr. energy repr. (— oo 5 oo)-"energy" repr. 

(162) 

where 

^ 2 = /res[S( z )] z=z „. (163) 

Physically, these links mean that the Gamow states yield a decay amplitude srf (zr — > £) 
given by the complex delta function, and that such decay amplitude can be approximated 
by the Breit-Wigner amplitude when we can ignore the lower bound of the energy, i.e., 
when the resonance is so far from the threshold that we can safely assume that the energy 
runs over the full real line: 

z n - E) = (-E\ Zn ) = iV^yKndiE - z n ) ^ . (164) 

V27T E-Zn 



Thus, the almost-Lorentzian peaks in cross sections are caused by intermediate, unstable 
particles. However, because there is actually a lower bound for the energy, the decay am- 
plitude is never exactly given by the Breit-Wigner amplitude. This means, in particular, 
that the standard Gamow states are different from the so-called "Gamow vectors" of [2]. 



Resonance expansions 

The scattering bras and kets are basis vectors that furnish the completeness rela- 
tion (4). The Gamow states are also basis vectors. The completeness relation (5) gener- 
ated by the Gamow states is called the resonance expansion. 

Resonance expansions are almost always obtained in the same way. One starts from 
the expansion in terms of bound and scattering states and then, by deforming the 
continuum integral into the complex plane, and by Cauchy's theorem, one extracts the 
contributions from the resonances that are hidden in the continuum and write them in 
the same way as the contributions from the bound states. 
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For the sake of simplicity, we will focus on the resonance expansion of the transition 
amplitude from an "in" state <p + into an "out" state <p : 

poo 

(<p-,<p+) = y dE((p-\E-)S(E)( + E\(p + ). (165) 

We now extract the resonance contributions out of (165) by deforming the contour of 
integration into the lower half plane of the second sheet of the Riemann surface, where 
the resonance poles are located, and by applying Cauchy's theorem. Assuming that the 
integrand ((p~\E~)S(E){ + E\(p + ) tends to zero in the infinite arc of the lower half plane 
of the second sheet, the resulting resonance expansion is 

00 p — 00 

(<p-,<p+) = £(<r|z„>(z„|<p + >+ / dE(<p-\E-)S(E)(+E\<p+). (166) 

n=\ Jo 

The integral in Eq. (166) is supposed to be done infinitesimaily below the negative real 
semiaxis of the second sheet. By omitting the wave functions in (166), we obtain the 
completeness relation (5). In Eq. (166), the infinite sum contains the contribution from 
the resonances, while the integral is the non-resonant background. Note that in (166) 
bound states do not appear, since the potential (78) doesn't bind any. 

In obtaining Eq. (166), we have assumed that the integrand (<p~|£ , ~)S(.E , )( + E , |<j!> + ) 
tends to zero in the infinite arc of the lower half plane of the second sheet. Since (166) 
makes physical sense, we are tempted to conclude that such must be the case. However, 
as showed in [9], ((p~\E~)S(E){ + E\(p + ) does not tend to zero in the infinite arc of the 
lower half plane of the second sheet. On the contrary, it diverges exponentially there. 
Therefore, Eq. (166) doesn't make sense as it stands. In order to make sense of it, one 
has to control the exponential blowup of (<p~\E~)S(E) ( + E\(p + ). The way to do so is 
by introducing an exponentially damping regulator e~ iaz , a>0. Thus, Eq. (166) should 
read 

00 p — 00 

(<p-,q>+) = \imYe- iaz »((p-\z n )(znW + )+ / dE e~ iaE (<p-\E-)S(E)(+E\(p+). 
o-o^i Jo 

(167) 

Physically, the regulator e Mz is simply the analytic continuation of the time evolution 
operator in the energy representation, e~' za = e~ lzt / h . Thus, the above regularized equa- 
tion must be understood in a time- asymmetric, time-dependent fashion as 

(9- = Y, e ~ iz4 ((p~\zn)(zn\(p + )+ / dEe- iEL *((p-\E-)S(E){ + E\(p + ) 

„=i Jo 

(168) 

for t > only. Equation (166) should then be seen as the (singular) limit of Eq. (168) 
when t — > + . That for resonances t = ah must be positive is in accord with the time 
asymmetry of (158). 

Expansions (167) and (168) are the reason why we chose a Gaussian falloff for the 
elements of 4> exp : When the wave functions have a Gaussian falloff in the position 
representation, we can regularize their blowup in the energy representation and interpret 
the regulator as a time- asymmetric evolution. 
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Resonance expansions allow us to understand the deviations from exponential decay. 
When a particular resonance, say resonance 1, is dominant, then the Gamow state of 
resonance 1 will carry the exponential decay, whereas the background, which includes 
in this case also the contribution from other possible resonances, carries the deviations 
from exponential decay. 

The full account of this lecture will appear in a forthcoming paper. 

CONCLUSIONS 

We have seen why the rigged Hilbert space, rather than the Hilbert space alone, is 
needed to formulate quantum mechanics when the observables have continuous and/or 
resonance spectra. The rigged Hilbert space captures the physics of continuous and 
resonance spectra better than the Hilbert space, because in the rigged Hilbert space 
physical quantities such as commutation relations, uncertainty principles and resonances 
have always a precise meaning. 

In addition to provide the mathematical support for Dirac's bra-ket formalism, for 
the Lippmann-Schwinger equation and for the Gamow states, the rigged Hilbert space 
can be used to obtain the resonance (decay) amplitude in terms of the complex delta 
function. Such decay amplitude can be approximated by the Breit-Wigner amplitude 
when the lower bound of the energy can be ignored. 

To finish, I would like to mention that there is still a long list of pending questions 
worth pursuing, such as the invariance properties of 4> exp under time evolution or a 
detailed proof of the asymmetry in the time evolution of the Gamow states. 
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